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Abstract. In the present work we introduce a novel multi-agent model with the aim to reproduce the 
dynamics of a double auction market at microscopic time scale through a faithful simulation of the matching 
mechanics in the limit order book. The agents follow a noise decision making process where their actions 
are related to a stochastic variable, the market sentiment, which we define as a mixture of public and 
private information. The model, despite making just few basic assumptions over the trading strategies of 
the agents, is able to reproduce several empirical features of the high-frequency dynamics of the market 
microstructure not only related to the price movements but also to the deposition of the orders in the 
book. 

PACS. Market Microstructure - Econophysics - Detrended Fluctuation Analysis - Multi- Agent Models 



2 1 Introduction 



In the past few years, following the increasing power of 
technological infrastructures, high-frequency trading, which 
broadly speaking includes every strategies which holding 
period is shorter than a day, has bloomed among the ma- 
jor financial institutions around the globe and, nowadays, 
it accounts for about 70% of all the volume traded in US 
equities. The same trend is evident also in Europe where 
the number of transactions for the most liquid contracts 
has recently experienced an exponential growth: in the 
Eurex Stoxx futures index, for example, the number of 
daily trades has gone from 12500 in Jan uary 2005 up to a 
maxim um of 150000 in November 2008 (|Bartolozzi et al.l . 
l2007aj) . As a consequence for the growing interest in the 
short time scales, the exchanges has started to provide live 
feeds of every single order submitted (buy or sell) , there- 
fore, making the market microstructure a primary field of 
interest for the financial practitioners. 

The progressive shift towards the very high-frequencies 
has not been unnoticed in the Econophysics community 
which, in the meanwhile, h as become an established field 
of study among physic i sts (iBouchaud and Potter;, 
Mantegna and Stanlev, 1999; ^Paul and Baschnage 
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1999; 



19991: 



Voit . 20051 . In the multi-agent framework, a common fea- 



ture in the physicist's approach to the market microstruc- 
ture of order driven markets is the lack of a proper utility 
function, which is often invoked in Economics/Econometric 
literatu re in order to study the rational behaviour of the 
agents (|Bailevl . l2005( ). and, consequently, a simple statisti- 
cal approach to the problem is preferred (JBak et al.L Il997l : 



telligence (JFarmer et al.L 120051) : a recent review on some 
of the models proposed so far can be found in ( Slaninal . 
2008). For an Econometric approach to the ma rket mi- 
crostructure, instead, the reader can refer to ( O'Haral . 
11997,) . 

Order driven markets are based on the principle of 
continuous double auction where the price of an asset is 
continuously adjusted in order to account for the inflow 
of demand and supply. In particular, the orders placed by 
the traders, that is the number of contracts they want to 
buy or sell along with the corresponding price, are orga- 
nized and matched in the limit order book (LOB) which 
ultimately defines the microstructure of the market. At 
every instant in time, the LOB, a sketch of which is given 
for clarity in Fig. [1] is described by two sets of limit or- 
ders in opposite "directions" , one for the long orders (buy 
side) and one for the short (sell side) . Each order is char- 
acterized by a limit price, that is the price the trader 
is willing to buy/sell, and the number of contracts re- 
quested/offered, or volume. The minimum price difference 
greater than zero between two orders is called tick size 
and, in real markets, depends on the specific contract. 
When more than one limit order is send to the same price 
then the exchange rank them in a stack by arrival time: 
the first to be executed will be the oldeso The mid-point 
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^ Note that while this ranking describes a typical situation, 
the execution of limit orders can slightly change depending on 
the rules of the specific exchange 



M. Bartolozzi: A Multi Agent Model for the Limit Order Book Dynamics 



of the simulations are presented in Sec. [3] Discussions and 
conclusion are left in the last section. 
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Fig. 1. Cartoon representation of the LOB. Note that gaps 
between price levels can be present, especially in non-liquid 
contracts. 



price, Pm, which is usually referred as the "price" of an 
asset, is the mean value between the best ask and the best 
bid, being, respectively, the lowest price on the ask side 
and the higher price on the buy of the LOB. The execu- 
tion of a limit order, that is the actual trade, is triggered 
if, and only if, an order in the opposite direction matches 
its price quote. The incoming order can either be another 
limit order, in which case the order executed is at best ask 
or bid, or a market order, that is a request to execute a 
certain volume, starting from the best price, immediately 
and until it has been completely filled. This latter type of 
order is considered very "aggressive" given that the price 
of the spread is paid upfront and it is usually associated 
with "impatience" traders. Cancellations also play an im- 
portant role in the dynamics of the LOB. In fact, limit 
orders that have not been filled in a reasonable amount 
of time, according to the trader's necessity, are usually 
removed from the LOB. 

In the present work we introduce a novel multi-agent 
modejj with the aim to reproduce the dynamics of the 
market microstructure. Our agents relate their actions to 
a stochastic variable, the "market sentiment", which in- 
cludes feedbacks from both public and private informa- 
tion. Besides, the model relies on a realistic LOB me- 
chanics which takes into account every possible change 
on "event" base such as the arrival of limit orders, market 
orders or cancelations as well as a faithful order matching 
mechanism. 

The paper is organize as following: in the next sec- 
tion we introduce our model while the numerical results 



For "agents" here we indicate all possible financial institu- 
tions such as investment banks, hedge funds etc... rather that 
the occasional trader. 



2 The Model 

2.1 Brief outline 

Our market model evolves in discrete time stepf|2| during 
which each agent may undertake a certain action or just 
wait for a more profitable opportunity. These actions can 
be broadly divided into cancellation and active trading, 
the latter including both limit and market orders. In par- 
ticular, all decision steps are based on dynamical proba- 
bilities which are function of both the private and public 
information, the former being related to the "state" of the 
market as we will see in the following sections. Moreover, 
we assume that our agents can have just one open posi- 
tion at the time and, therefore, before issuing a new order 
they need to close their current position. By using this 
limitation we neglect explicit market making activity. 

Another central part of our model is the order gener- 
ation step where the specifications for each order such as 
type (limit or market), price (for limit orders) and volume 
are decided. In the next sections we explain in details all 
the fundamental building blocks. 

2.2 Public and private information update 

At the beginning of each time step, we assume that the 
agents have access to the current state of the LOB, that 
is they can "see" every order placed in the market: all 
the indicators derived by this knowledge, such as the mid- 
point price for example, are classified as public informa- 
tion. However, these are not used in their "raw" form but 
they get "smoothed" : this filtered information, that we re- 
fer as "perceived" , is the one which is actually used in the 
process of decision making. The reason behind this extra 
step is to take our simulations one step closer to reality 
where indicators are usually pre-processed in order to get 
rid of some noisqj. The filter used in our simulations is an 
exponential moving average (EMA) (JFusai and Ronconil . 
2008^ ) which is defined for a generic time series x{t) as 



m = ^x{t) + fi-i]x{t-i) 



(1) 



where L is proportional to the memory of the process. 
By using the previous equation Eq.((T]), we define the per- 
ceived volatility as 



m = ^/r^t), 



(2) 



It is important to stress that in real markets the arrival 
of new information at microscopic time scale is heterogeneous 
in time and often people consider the trade time as a more 
appropriate scale. While the output of the model is updated in 
discrete time steps for simplicity, the LOB is changed on event 
base. 

* As a consequence, short term memory is induced in the 
signal itself. 
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being r the one step return of the mid point price 
r{t) = Pm{t)—Pm{t—1)- Note that, for time being, we have 
not included in the public information any news realize, 
the i mportance of which has been question ed over time , 
see (jCutler et all . Il989[ ) for example or poulin et all . 
I2OO8I) for a more recent criticism. 

Regarding the private information, instead, we assume 
that it can be represented by a simple Gaussian process, 
independent for each trader, with zero mean and standard 
deviation proportional to the perceived volatility of the 
market, vlt): this information can be thought to represent 
the convolution of the trading strategies used by the agent. 



2.3 Cancellation of orders form the LOB 

At each time step, agents having an outstanding order 
in the LOB will evaluate the possibility of removing it 
according to two different criteria. The first corresponds 
to a simple time-out, that is the order is automatically 
removed if it has not been executed in a Tmax number 
of time steps relatively long if compared to the average 
time for a transaction and fixed to 100 time steps in our 
simulations. The second criterion, instead, is related to a 
strategic decision based on the current market condition. 
In particular, we define a cancellation probablity, ^Jc,{t) G 
[0, 1], common to all the agents as 



V'^W-1 



-yi>{t) 



(3) 



being 7 = 0.02 a sensitivity parameter. The former ex- 
pression, which relies just on the perceived volatility h'{t) 
in terms of information, is justified by the empirical ob- 
servation that the cancellation frequency increases when 
the market is highly volatile: this is clearly a case of self- 
reinforcing effect tri ggered by the herding behaviour of the 
market participants (jCont and Bouchaudll2000tlBartolozzi 
12004 . 



2.4 Active trading and market sentiment 

While the cancellation step takes place, agents with no 
orders in the LOB evaluate the possibility to enter the 
market. Their decision is based on a stochastic variable 
which represents the "level of confidence" in their price 
forecast: the market sentiment, cf). This quantity, which 
is a core part of our model, relates the public and the 
private information through a multiplicative process and, 
specifically, for the zth agent at time t we define 



— The parameter Ki represents the agent's degree of risk 
aversion and it used in order to diversify the appetite 
for risk. Its value is randomly selected, at the beginning 
of the simulation, from a uniform distribution bounded 
in [0.25, 0.75]. 

— Tpp — l — ^v{i) is the volatility risk and mimics the fact 
that traders are more cautious to send orders during 
high volatility periods given that the risk associated to 
these is higher as well. 

— i]{t), instead, is a proxy for the liquidity risk: a sparse 
LOB has a potential large negative impact on the exe- 
cution of a trade and, therefore, it decreases the prob- 
ability that an agent is willing to accept the risk. In 
the present work we assume r]{t) = N{t)/N where A'^ is 
the maximum number of orders in the LOB, equivalent 
also to the total number of agents in the simulation, 
and N{t) the number of orders in the LOB at time 
step t. From the definition we have that ri{t) G [0, 1]. 

— Ei {t) represents the private information which is drawn 
from a Gaussian distribution with zero mean and stan- 
dard deviation equal to 0[t), ei{t) = i>{t) ■ G{Q, 1). 

The market sentiment, Eq.Q, can be thought as the 
convolution between the agent's trading strategies, the 
private information, and the risk factors evaluated via the 
public information (volatility and liquidity in this case): 
the stronger is the signal the more likely will be for the 
trader to take a decision. 

The next step involves the mapping of 0^ [t) into a trad- 
ing probability, Pi(t) G [0, 1], via a the following transfer 
function 

p,(t) = -|arctg[0,(i)]|, (5) 

TT 

which represents the probability for the ith agent to sub- 
mit an order at time step t. Conversely, an agent will not 
take any action with probability 1 — Pi(t). 

The dire ction of a trade, assuming that short selling 
aiidaHaaTOad llowed. di{t) = -1-1 for long orders (buy) and 
di{t) = — 1 for short orders (seU), is also derived from the 
market sentiment according to 



d,it) 



.it) 



\MtW 



(6) 



(/),(i) =(/)o -K, -ipiit) -riit) •e,(i). 



(4) 



where each terms represents a different aspect which may 
impact on the decision making process as following: 

— 00 is a strength parameter common to all the agents 
which, as we will see in the next section, fixes the trad- 
ing frequency of the model and, therefore, the time 
scale. 



as if 0i(i) was a momentum indicator for the market 
direction. Besides it is worth underlying that, according 
to the previous definition, the direction of a trade depends 
only on the sign of the private information, Eq.Q, and, 
therefore, there is 50% chances for a trade to be a buy or 
a selE 

2.5 Order generation 

Each time an agent submits an order to the market its 
specifications are defined in order generation step which 
addresses the following two points: 

^ The market sentiment, (j>i{t), and consequently the trade 
direction, can be statistically skewed in one direction depend- 
ing on different market factors and, therefore, leading to herd- 
ing effects. The study of this phenomenon will be the topic of 
future investigations. 
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— Limit or market order? The decision to submit a limit 
or a market order is related to the impatience of the 
trader: if an order needs to be filled as soon as possible 
he/she may accept to pay the cost of the spread up- 
front and send a market one. The alternative would be 
a less aggressive limit order that will execute just at 
a pre-fixed price: the further away from the opposite 
best the more time the order will require to get filled. 

— Size of the order? The number of lots that the trader 
is willing to buy or sell can be related to several fac- 
tors such as to the specific execution strategy, to the 
available liquidity in the market, the volatility at that 
point in time etc... Moreover, in a real trading envi- 
ronment single large orders are usually splitted into 
smaller ones in order to minimize the cost related to 
their impact. While the latter is an important issue 
in practice, in order to keep things simple, we do not 
tackle this problem in the present work: orders are sent 
in just one single chunk. 

In our model the first point is addressed by setting 
the submission price of an order, which can be interpreted 
as the degree of "aggressiveness", probabilistically via a 
value drawn from log-normal distributior{3, ^, as following 

Pi^Pb-i^~ Qp), (7) 

for a long order (buy) Pi while for a short (sell), P, 



LinMOider Market Order 




Price 



Fig. 2. The price at which an order is submitted, and therefore 
its "aggressiveness", is drown from a log- normal distribution 
which g-quantile is centered at the best bid (long order) or best 
ask (short order). The "direction" of the distribution is also 
chosen in accordance with the order direction. In the cartoon 
above we sketch the probability of price submission for a long 
(buy) order. 



Ps^Pa + {^-Qp), 



(8) 



being Pb and Pa the best bid and the best ask, respec- 
tively, and Qp the g-quantile, q = 0.5 in our case, of the 
distribution. Besides, the value of ^ is rounded up to an 
integer value and, therefore, fixing, without losing of gen- 
erality, the tick size of our market to 1. It is also important 
to stress that following the former procedure we explicitly 
assume that limit orders can be submitted asymptotically 
far from the best price as shown in the sketched in Fig. [51 

After the submission price has been fixed, the type of 
the order is decided based on its relative position to the 
best prices: if the submission price results to be greater 
than the ask price and the order is long (or lower than 
the bid price and the trade is short) then we interpret this 
as a market order and all the volume will be completely 
filled starting from the best price. All the other orders are 
considered limit orders and, for a fixed price, they are or- 
ganized in stacks from the oldest to the newest being the 
formers the first to be executeqj. This modelling approach 
is motivated by two empirical observations: the distribu- 
tion of the volumes in the LOB, when averaged over long 
periods, displays "fat tails" and that the frequency of mar- 
ket order represents just a fraction of that ofjimit orders 
submitted in the market (JBouchaud et al.L l2002l) . 



The parameters used in the simulations for the log-normal 
distribution are 7 for the mean and 10 for the standard devia- 
tion. 

"^ If the submission price coincides with the opposite best 
we give 50% probability for the trade to be a limit order or a 
market order. 



The second issue concerns the order size, that is the 
number of contracts, or volume, that an agent is willing 
to buy or sell. Coherently with the stochastic nature of 
this process, this is drawn from a log-normal distributior{3, 
rounded to integer, with the following constrains: the or- 
der has to be greater or equal than one and smaller than 
one quarter of the total volume contained in the appro- 
priate side of the LOEJj. 

Before moving to the next section, we wish to under- 
line that the order generation step used in our model is 
just a first order approximation of what really happens 
in real execution strategies. In particular, we do not ad- 
dress the problem of splitting large volumes into smaller 
ones, typically used in order to minimize the transaction 
costs: this issue, despite being very important for market 
impact, goes beyond the scope of this work at present. 



The parameters of the distribution used in the simulations 
are the same as for the price submission that is 7 for the mean 
and 10 for the standard deviation. 

^ In the simulation this limit is hardly reached: usually, single 
orders tend to be relatively small compared to the total volume 
available in the LOB. Moreover, before submitting a market 
order we also perform a liquidity check: if the number of lots 
in the requested side of the LOB is less or equal Nmin, with 
Nmin = min(50, N/10), then the order is not issued. Liquidity 
checks on the volumes in the book are often done also in real 
trading in order to have an estimate of the risk associated with 
the trade itself. 
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Fig. 3. Kurtosis of the one step returns estimated over an 
ensemble of three simulations each one consisting in 10^ time 
steps. It is noticeable how this quantity converges to 3, the 
value for a Gaussian process, as (po increases. By decreasing 
00, instead, the activity gets very low and the dynamics of 
the model starts to be ruled by few large market movements. 
The parameters for these simulation have been declared in the 
previous section except for the EM A history L, Eq.([l]), which 
in this instance has been fixed to 5. 
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Fig. 4. (Top) A window of the one step returns time series 
generated by the model with N = 10000, ^o = 0.165 and 
L = 5. The time series, standardized according to r{t) ^■ 
(r(t) — f)/a{r) being f the average return and a{r) the stan- 
dard deviation, displays an intermittent character. (Bottom) 
Autocorrelation function, p{t), related to the time series on 
top: a significant anticorrelation is evident for the first time 
steps. 



3 Numerical simulations at short time scales 

In this section we report the results of the numerical sim- 
ulations of our microscopic market model for N — 10000 
agentq^. The trading frequency, which can be interpreted 
as a time scale for the simulation, is fixed by the param- 
eter (/)()■. as its value gets smaller the trading activity be- 
comes relatively lower and gaps within the LOB, leading 
to large price changes, are more likely to appear. On the 
other hand, if the activity is very high, the book is almost 
always full and large fluctuations will occur more sporad- 
ically. The scaling of the fluctuations with the activity is 
emphasized in Fig. [3] where we report the kurtosis, k, of 
the one step returns against (/jq. 

For 00 — 0.165, the value that we use in the rest of 
this section, the activity resembles that of the market at 
very short time scales, from seconds to minutes depending 
on the specific contract, where non-Gaussian fluctuations 
play a fundamental role. In this regime, which is the one 
of interest in the current work, there is a probability of 
approximately 30% for the mid-point price to remain un- 
changed after one time step. 

Another important parameter is the EMA history, L 
in Eq.([T]). In fact, we have found that activity clustering, 
such as high volatile periods, are particular evident when 



^^ Numerical tests have shown that the outcomes are indiffer- 
ent to the number of agents as long as their number is relatively 
large. 



3 < L < 10. In the simulations we have fixed L = 5 for all 
the agents. 



3.1 Price and volatility dynamics 

A sample of the time series of one step returns, r{t), 
generated by the model is reported in Fig. |4] along with 
its autocorrelation function, p(r). From the plots we can 
notice an intermittent dynamics, characteristic of finan- 
cial time series at short time scales, as well as a signifi- 
cant negative correlation up to few time steps. This lat- 
ter effect, known empirically as "bid-ask bounce", in the 
present model is purely a result of the order book me- 
chanics. However, for real high-frequency data this effect 
can l ast up to few minutes depend ing on the specific mar- 
ket (jBouchaud and Potterd . I1999D : this is an important 
indication that, on top of the LOB mechanics, there must 
be other mechanisms, such as memory feedbacks or order 
splitting for example, responsible for the enhancement of 
the negative correlation in price changes. 

The anticorrelated behaviour of the returns is also con- 
firmed by an estimate of the Hurst exponent, H, done via 
the detre nded fluctuat i on ana lysis algorithm ( Peng et al.l . 
.1994; Ba rtolozzi et al.l . [2007b() . According to this method, 
origina l ly developed durin g the Nile's river dam project 
(JHurstl . I1951I : iFedeii Il988l) . a time series is persistent if 
H > 0.5, anti-persistent ii H < 0.5 or uncorrelated ii H = 
0.5. The value found for the one step returns in our model 



M. Bartolozzi: A Multi Agent Model for the Limit Order Book Dynamics 



10"' 


A 


- 


10-'' 


/ \ 


: 


^ 10"' 


: \ 


: 


10"* 


j \ 


: 


10-' 


/ \ . 


♦ : 




'''r*Wi)', 



' ' m, ,:» . 



I' ,,, ',1 



-15 -10 



10 15 



Fig. 5. Probability distribution function, !?■, for the one step 
returns standardized as in Fig. [4] (solid line). This distribution 
is leptokurtic, that is the tails are "fatter" when compared to 
those of a Gaussian (dotted line) . 



]s H — 0.495(2) over an ensemble of 3 runs of 10^ sam- 



Fig. 6. (Top) Absolute value of returns, or instantaneous 
volatility, for TV = 10000, (t>o = 0.165 and L = 5. Patches of 
high-volatility periods are visible. (Bottom) Autocorrelation 
function relative to the time series on top. The instantaneous 
volatility results to be correlated positively up to few hundreds 
steps. 



pies each, being the error on the last digit, re ported in the As p ointed out in the previous section, we would like 

brack ets, estimated via the bootstarp method ( Efron and TibbhialaBii s the role played by the memory induced via the 



[l99J). Noticeably, by removing the zero from the time se- 
ries the value of the former exponent is statistically equiv- 
alent, in fact we find in this case Hq = 0.491(3) 

The probability distribution function or pdf^ denoted 
as !?", of the one step returns is reported in Fig. [S] where 
the large fluctuations observed in Fig. |4] (top) give rise 
to a leptokurtic shape of the distribution, that is the tails 
are "fatter" than those of a Gaussian. This feature, per- 
sistent up to time scales o f weeks, is well documented 
in different em pirical works ("Bouchaud and Pottersl . 119991 : 
iMantegna an d Stanley, 1999; Voit, 2005|). However, while 
in most of the previous examples the asymptotic decay of 
the tails can be described by a power law, in our simula- 
tions this is more consistent with an exponential one. From 
the same plot it is also possible to notice a relatively large 
presence of zero returns, which is also a common feature 
of financial price time series at very high-frequencies. 

In Fig. [5] (Top), instead, we report the instantaneous 
volatility defined as u{t) = \^{t)\ where it is noticeable 
the presence of clustering, that is periods of high activ- 
ity tend to be patched together. This memory effect is 
highlighted by the slowly decaying autocorrelation func- 
tion. Fig. [5] (Bottom), and by the Hurst exponent estimate 
which gives the value oi H = 0.61(2), indicating a signifi- 
cant persistency. Besides, we have verified that the corre- 
lation, despite being weaker, is still present after the zero 
returns removal, Hq — 0.56(1). A similar behav iour is also 



perceived volatility, Eq.([5]), and, therefore, the parameter 
L, in the volatility clustering. In fact, the clustering phe- 
nomenon seems to be particulary relevant only when the 
former feedback is present: this is a strong indication that 
the use of short moving averages in high-frequency trading 
strategies, implicitly inducing memory effects, can be the 
source of the bursts of patched volatility observed at short 
time scales. Also note that, following the previous observa- 
tion, in our model both the market sentiment, Eq.(j31), and 
the cancellation process, Eq.([2]), contribute to the volatil- 
ity clustering. In fact, it can be shown via numerical ex- 
periments that by "switching off' the dependence from 
i>{t) in one of these two terms, for example by substitut- 
ing it with a constant, the clustering effect gets noticeably 
reduced. 



3.2 Traded volume and impact 

So far we have examined the bahaviour of the price returns 
and volatility. However, both these quantities are related 
to other fundamental factors, one of the most important 
being the traded volume, V. Some authors, in fact, have 
recently argued that large transactions could be responsi- 
ble for the non-Gaussian fluctuatio ns observed in t he fat 
tails of the returns distribution (Ga baix et al.l . |2003| ) even 
though more recent empirical work seems to suggest that 



liquidity crises, that is the temporary pr esence of gaps in - 
obser ved in real markets, see, for example , in (iLiu et al.l . the L OB, should be at the origin of those (iBotichaud et al. , 

Il997t iGoDikrishnan et all Il999l : iLhTet al.l . I1999Il 120021: iDaniels et al.l . 120031: iFarmer et all . 120041: iLillo and Farmed . 
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Fig. 7. (Top) Window of 3000 samples for the traded volume Fig. 8. Instantaneous impact function where the negative part 
in the model along with the relative autocorrelation function corresponds to sell initialized orders. Note that while the im- 
(Bottom). The parameters used are A*' = 10000, (jio = 0.165 pact is linear for V ^ 0, for large volumes this feature is lost, 
and L — 5. Also, the values of the traded volume in the plot have been 

rescaled by the maximum volume traded in order to fit the in- 
terval [-1 1]. In this case the rescaling has been made simply to 
20051: IWeber and RosenowL [20051 : lFarineil . l2006HWeber and Ifegifetento the visual comparison of different impact functions 



20061) . Moreover, the dynamics of the traded volume time 
series, like the volatility, has been observed to be charac- 
terized by clus ters of intense activity fo llowed by relatively 
quite periods (JLillo and Farmen . l2004l) . 

For our model, the cumulative traded volume during 
each time step is reported in Fig.[7J where a behaviour very 
similar to that of volatility is evident. The intuitive rela- 
tion between the two is further confirmed by their cross- 
correlation coefficient which we found greater than 30%. 
Moreover, the estimate Hurst exponent for the traded vol- 
ume results H = 0.70(2) (iJg = 0.66(1)), matchings rea- 
sonab ly well the value found empirically in (jLillo and Farme: 
120041 ■ namely H = 0.73(7). It also worth pointing out that 
a possible introduction of order splitting in the execution 
part of the model, which is missing at the moment, should 
enhance this correlation. 

We also estimate the average instantaneous change in 
price subsequent to a trade of a certain volume V^ usually 
known as impact function. This quantity is very important 
in practical applications being a proxy for the liquidity of 
the market: its value indicates how, approximately, an or- 
der can penetrate "deep" into the LOB. The impact func- 
tion for our model. Fig. [8j displays a linear behaviour as 
long as V is relatively small. However, due to the presence 
of "volume barriers" , as we will see in the next section, 
this proportionality is lost as the volume increases. While 
the same sha pe for the impact function is obser ved in real 
markets, see (jWeber and Rosenowl . l2005l l2006l) for exam- 
ple, zero-intelligence models usua lly reproduce only the 
linear part, with the exception of ( Farmer et al.l . l2Q04[ ). 

The values reported in Fig. [5] have also another impor- 
tant implication. In fact, the average price response to the 



in case of changing the parameters of the volume distribution, 
see Sec. 12.51 The simulation has been run with N — 10000, 
4)0 = 0.165 and L — 5. 



large volumes is not enough to justify the extreme fluctu- 
ations observed in the pdf of returns that can be greater 
than 10 ticks. Fig. \S\ which, therefore, have to be related 
or to rallies of orders in one direction or to a temporary 
lack of liquidity between the levels of the LOB. Of course 
nothing prevents both scenarios to be realized at the same 
-itime. 

a 



3.3 LOB average shape, spread and imbalance 
dynamics 

In this final section we turn our attention to the way the 
orders are "on average" deployed in the LOB as well as 
the dynamics of some quantities related to the relative 
position of the orders such as the spread and the imbalance 
between buy and sell. 

Firstly, we show the "mean" shape of the LOB, that is 
the average volume present at a fixed distance, \P — Pm\, 
from the mid-point price. From the plot. Fig. [51 it is clear 
that the liquidity increases steeply up to a maximum, lo- 
cated few ticks away from P,„, and then it starts to de- 
cay as we move away from it. The peak in volume can 
be thought as a sort of "volume barrier" which prevents 
large orders to get too deep into the LOB and, therefore, 
confirming the speculations made in the previous section. 
Very similar shapes have been found in empirical studies 
where it has been suggested that the asymptotic shape of 
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Fig. 9. Average shape of the limit order book reported against 
the distance from the mid-point price. Note that, given the 
symmetry of the model, we have aggregated the data on the 
ask and the bid side of the book so to have better statistics. 
The parameters used are the same as in Fig. [T] 



Fig. 10. (Top) A sample for the spread time series (in ticks): 
it is possible to notice how, occasionally, the spread can be- 
come relatively large. (Bottom) Autocorrelation function for 
the spread time series. The parameters used are the same as 
in Fig. El 



the LOB, related to the "patient " trad ers, may follow a 
power law decay (JBouchaud et al.l . l2002f ). It is also impor- 
tant to stress that the smooth "mean" shape of the LOB 
is not significant for the "instantaneous" shape which, in 
fact, can be relatively sparse. 

Another important quantity related to the placement 
of the orders in the LOB is the spread, S. This is the 
difference in price between the best ask and the best bid 
and it represents the cost that a trader has to pay up- 
front in order to execute a market order. Moreover, the 
spread is a further proxy for the liquidity: large spreads 
would indicate shallow markets while very liquid ones will 
tend to keep a very small spread at all times, possibly 
close to one tick. Dynamically, the sp read displays persis- 
tency in time with 7 < H < 0.8 (iPlerou et al.l 120051: 
ICaiueiro and Tabakl 120071: iGuet all . 120071 ) depending on 
the market and the time scale of observation, while the 
asymptotic shape of its pdf can be ch aracterized by a 
power law function (jPlerou et al.L l20Q5ll . The same two 
features appear also in our simulations, Fig.[TUl even though 
the Hurst exponent result to be relatively smaller if com- 
pared to the empirical findings, namely H = 0.56(2) {Hq — 
0.56(2)): this discrepancy may arise from the fact that at 
the moment we are considering only trivial execution al- 
gorithms. 

Lastly, we focus our attention to the dynamic imbal- 
ance between buy and sell orders or simply the volume 
imbalance, AV{t), defined as 



rb.a 



AV{t) 



Nbit) 






(9) 



i=l 



i=l 



where, for a time t, V^ ' represents the volume of ith limit 
order and Nf,^a their total number, respectively, for the bid 
and the ask side of the LOB. The time series, Fig. [TTJ as 
expected for a market near-equilibrium, displays a mean 
reverting behaviour, that is the imbalance tends to fiuctu- 
ate all time, more or less symmetrically, aroun d a "pseudo- 
equilib rium" price, possibly close to criticalitv'Bartoloz zi et all 
(j2005i |2006() . The same dynami cs has been observed for 
AV in different futures contracts ([Bartolozzi et al.l . l2007af) . 



4 Discussion and conclusions 

In the present work we have developed a multi-agent frame- 
work characterized by a realistic order book keeping as 
a tool for the study of the activity of a double auction 
market at microscopic time scales. Inside this framework 
the model relies just on few basic assumptions related to 
the agent's strategic behaviour. Among the most impor- 
tant ones, the order submission process makes use of a 
stochastic variable, the market sentiment, which is related 
to both the public and the private information. Besides, 
the agents lack of the concept of "utility maximization" 
often invoked in the economics literature and embodying, 
in mathematical terms, the concept of "rationality" . 

Despite its simplicity, the model manages to reproduce 
several empirical features of the high-frequency dynamics 
of the stock market such as the negative correlation in 
market returns, clustering in the trading activity (such as 
volatility, traded volume and bid-ask spread) as well as 
a non-linear response of the price change to the volume 
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Fig. 11. (Top) Time series of imbalance between demand and 
supply displaying a clear mean reverting nature. The parame- 
ters used for the simulation are the same as in Fig. [T] 

traded. Moreover, the similarities with the real markets 
extend also in the way the orders are deployed on the 
LOB as observed through the average shape of the book 
and the volume imbalance time series. 

In conclusion, our model indicates that large part of 
the dynamics of the stock market at very short time scales 
can be explained without the requiring any particular ra- 
tional approach from an agent prospective if not some 
memory feedback which, in our model, are represented 
by short term moving averages of the public information. 
Moreover, our results confirm that large price movements 
are more likely to be related to a temporary lack of liquid- 
ity in the LOB rather than to large volume transactions. 

Our future work will involve adding more realistic fea- 
tures and feedbacks in the model. In particular, one in- 
teresting effect to take into consideration would be the 
herding phenomenon which, it has been argued, to be at 
the origin of dramatic liquidity crises which can ultimately 
lead to financial crashes (jSornettd . |2003). Other factors 
that can claim to play a fundamental role at these time 
scales are news realize and market making strategies: both 
of them should to be taken into account. 
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